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Abstract 

Coherent charge transfer in a hnear array of tunnel-coupled quantum dots, 
electrostatically coupled to external gates, is investigated using the Bethe 
ansatz for a symmetrically biased Hubbard chain. Charge polarization in this 
correlated system is shown to proceed via two distinct processes: formation of 
bound states in the metallic phase, and charge transfer processes correspond- 
ing to a superposition of antibound states at opposite ends of the chain in 
the Mott-insulating phase. The polarizability in the insulating phase of the 
chain exhibits a universal scaling behavior, while the polarization charge in 

the metallic phase of the model is shown to be quantized in units of e/2. 
PACS numbers: 73.20.Dx, 73.23.Hk, 71.10.Fd 
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Tunneling of a single electron from one region to another in a mesoscopic system leads to 
a modification of the dielectric response of the system [|l| which can be detected via single- 
electron capacitance spectroscopy 0. Capacitance measurements allow one to study charge 
transfer in equilibrium, and thus provide an important alternative to transport measurements 
in probing the effects of coherent tunneling. In this Letter, we investigate the charge 
polarization of a linear array of tunnel-coupled quantum dots embedded between the plates 
of a capacitor (Fig. 1). The quantum corrections to the classical two-terminal capacitance 
of the system are shown to exhibit sharp resonances whose structure reveals directly the 
spatial correlations of the interacting many-body ground state of the system. We find 
that the localized character of the many-body states in the Mott-insulating phase of the 
model leads to extremely sharp capacitance resonances, which obey a universal scaling form 
analogous to the conductivity of the system On the other hand, the extended quantum 
states in the metallic phase of the model are shown to lead to fractional charge transfer, in 
contrast to the integer charge transfer predicted in Ref. 

The integrability of ID quantum many-body systems with open boundary conditions 
was first established p] for the one-dimensional (ID) Hubbard model. Bethe ansatz solution 
was recently extended to include boundary potentials , and the spectrum of bound states 
for a single attractive boundary potential has very recently been investigated |]IU[. Here, we 



investigate a Hubbard chain with equal and opposite boundary potentials at each end, which 
serves as a model of a capacitively biased ID array of quantum dots. In addition to the 
bound states found for the case of a single boundary potential [|10[, we find charge transfer 



states, which are quantum mechanical superpositions of antibound states at opposite ends 
of the chain. These charge transfer states are shown to dominate the polarizability in the 
Mott-insulating phase of the model. 

We consider a closed linear system of quantum dots coupled electrostatically to bias 
gates and a backgate (Fig. 1). The backgate allows the system to be charged with N 
excess electrons, this excess charge being shared among the dots in the chain by quantum- 
mechanical tunneling. We describe this coupled quantum dot chain by the Hubbard model 
|Tl| in the experimentally accessible limit when the interdot capacitances are negligible 
compared to the capacitances Cg to the external gates. The Hamiltonian of the system, 
including the work done by the external voltage sources, is 

H = -tY: E (cli.Q. + H.c.) + ^ E p,^ - CoVy2 + ^ (PL - Pi) , (1) 

cr i=l ^ i=l ^ 

where creates an electron of spin a on dot i, pi = Y^acl^Cifj, U = /Cg is the charging 
energy of a quantum dot, and Cq is the classical geometrical capacitance between the left 
and right gates. Eq. (1) is the prototypical minimal model of correlated fermions on a 
lattice, and describes e.g., the correlation-induced metal- insulator transition 0. The new 
feature investigated here is the nonperturbative effect of the external bias {V) described 
by the last term in Eq. (1), which polarizes the system. Unlike previous investigations of 
the charge response of the system [0, we do not treat the bias V as a. weak perturbation, 
but consider arbitrarily large values of V , leading to a finite transfer of charge across the 
chain. The polarization charge Q induced on the external capacitor plates characterizes the 
measurable dielectric response of the system. At zero temperature, the expectation value of 
the polarization charge is given by 
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{Q) = {QL-QR)/2 = -dEo/dV, 



(2) 



where Ql (Qr) is the polarization charge on the left (right) capacitor plate and Eq is the 
minimum eigenvalue of Eq. (1). The two-terminal capacitance of the device is defined as 
Cfj, = —d'^Eo/dV'^. These quantities can be exactly obtained for the quantum dot chain 
using the Bethe ansatz technique, as described below. 
The eigenvalues of Eq. (1) may be expressed as 

N 

E = -2tJ2 cos kj - CoVy2, (3) 
i=i 

where the pseudomomenta kj are a set of distinct numbers which satisfy the coupled 
equations 

c N i2k,{L+i) ^ TT sinkj- Xp + iU/At sinkj + + iU/At 
^ ^ ^ ^ ji^ sin kj -Xp- iU/U sin kj + Xp - lU/At ' ^ ' 



^K- sin kj + iU/At Xa + sin kj + iU/At _ Xa - Xp + iU/2t + + iU/2t 

/=i - sin kj - iU/At A« + sin kj - ill /At ~ Xa - Xp - iU/2t A„ + A/3 - iU/2t ' 

(5) 

where Aq, a = 1, . . . , M are a set of distinct numbers referred to as spin rapidities, and 



1 - {eV/2t)^e~'^'^ 



is the single-electron scattering matrix associated with the boundary potentials. 

The capacitive response of a chain of 4 quantum dots, calculated from Eqs. (2-6), is 
shown in Fig. 1 for several values of A^. The polarization induced by the external bias 
V leads to a transfer of charge across the system, which is reflected in the appearance 
of complex roots of the Bethe ansatz equations, corresponding to bound and antibound 
states on the boundary dots (see Table I). Let us first consider the Mott-insulating phase 
of the system, which occurs |l7|JTT[| for commensurate electron density, N = L. For low 
bias, eV < 2t, the Bethe ansatz ground state contains only real pseudomomenta, and the 
charge distribution remains essentially symmetric. For 2t < eV ^ U, a. bound state forms 
on the leftmost dot, characterized by the complex pseudomomentum kL- However, due to 
the incompressibility of the Mott insulator, an antibound state on the rightmost dot is also 
filled (fci_i), and there is thus no net transfer of charge. The Mott-Hubbard gap is reflected 
in the suppression of the low-bias capacitance (the dash-triple-dotted curve in Fig. 1). For 
a bias larger than the Mott-Hubbard gap, however, it becomes energetically favorable to 
depopulate the antibound state on the rightmost dot and populate an antibound state in 
the upper Hubbard band on the leftmost dot (region II in Table I). The pseudomomentum 
of this antibound state contributes — 2tcosA;L_i = [{U — eV/2)'^ + At"^]^^"^ ^ f/ — eV/2 to 
the ground state energy in Eq. (3) (plus small backfiow terms), indicating the presence 
of a second electron on the leftmost dot. The resulting transfer of an electron across the 
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array leads to a sharp capacitance resonance at eV = eVi ~ f/ in Fig. 1. Finally, for 
eV > eV2 — 2U, this antibound state becomes a bound state. 

In order to elucidate the nature of the charge-transfer resonance in the Mott insulator, 
let us first consider the simplest case L = 2, for which Eq. (1) reduces to a simple 4x4 
matrix. The polarization charge and capacitance may then be obtained directly [neglecting 
terms of order (t/f/)^], 

Q-CqV ^ 1 ^ 1 eV-U 

e 2 2^st^^{U-eVf 



Eqs. (7) and (8) predict a charge transfer of e across the chain and a capacitance peak at 
eV = U. Eq. (8) was obtained previously in Ref. where it was shown to describe charge 
transfer between two arbitrary mesoscopic systems coupled weakly by tunneling. Following 
the above argument on the nature of charge transfer in the Mott insulator, one may expect 
a result analogous to Eq. (8) to hold for larger chains as well, since the effective coupling 
of the boundary dots via the intervening Mott-insulator should decrease exponentially with 
system size. Indeed, the capacitance peaks at eV ~ U are found to become increasingly high 
and narrow (the area, which corresponds to the total charge transferred, is conserved) as L 
increases, but their shape is found to be described very well by Eq. (8), with t replaced by 
an effective charge transfer matrix element tes, as shown in Fig. 2(a). Fitting the calculated 
capacitance to Eq. (8), the effective charge transfer matrix element is found to have the form 

as shown in Fig. 2(b), where the correlation length ^ in the Mott-insulating phase of the ID 
Hubbard model is given by [0 



4t f°° \n{y + 



l/aU/t) = - dy '" 'I ' (10) 

Eq. (9) indicates that the effective charge transfer matrix element, which characterizes the 
resonant polarizability of the Mott insulator, exhibits a finite-size scaling analogous to the 
conductivity of the system, which also decreases exponentially with system size in the Mott 
insulator 0. tefr is in fact related to the equal-time Green's function, tes = tG{l,L) = 
tJ2a{^\ciaCLa\0) , and it has already been argued that G has the same scaling form as 
the conductivity for another choice of boundary conditions. Dielectric measurements thus 
present the intriguing possibility to study experimentally the correlation length of a Mott 
insulator formed in a coherent system of quantum dots. 

While the charge transfer resonances in the Mott-insulating phase of the model are well- 
described by the theory of Ref. it is evident from Fig. 1 that the capacitance in the 
metallic phase of the model, which may exhibit a low-bias double peak structure, can not in 
general be described by an equation of the form of Eq. (8). As shown in Fig. 3, this double 
peak structure in the dilute metallic phase of the system is accentuated with increasing 
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system size, and corresponds to a polarization charge with well-defined plateaus quantized 
in units of e/2, unlike the integer charge transfer described by Eq. (7). From Table I, we see 
that the polarization of the system in the metallic phase proceeds via the successive capture 
by the boundary dot of electrons from the Luttinger liquid states delocalized along the chain 
(since the antibound states are empty for N < L), the first at eV = 2t and the second at 
eV = eV2 — 2U. The breakdown of Eqs. (7) and (8) is due to the fact that the system can 
no longer be divided into just two weakly coupled subsystems, as was assumed in Ref. [|I|, 
and instead becomes one coherent whole in the metallic phase. The fractional increments 
of polarization charge shown in Fig. 3 arise because the trapping of an electron from the 
gapless quantum states in the central part of the array leads to an effective charge transfer 
over only half the system. As a result, the magnitudes of the charges induced on the 1-st 
and L-th dot are not equal to each other in the metallic phase in contrast to the equal and 
opposite charge polarization at both ends of a quantum dot chain in the insulating phase. 
For example, for = L — 1 electrons in the chain, the polarization charge on the L-th dot 
maximizes to one unit of charge — e at eV ~ 2t, whereas the polarization charge on the 
1-st dot increases by +e at eV ~ 2U. The plateau structure in the capacitance near zero 
bias as seen in Fig. 3 arises in the dilute quasiparticle/quasihole metallic limit (A^ ~ 1 or 
L — N ^ 1), and its disappearance with increasing filling is an experimentally observable 
strong correlation effect in the charge polarization of the metallic phase. 

Let us comment on some of the idealizations employed in the above calculation. The 
introduction of an interdot capacitance, neglected in Eq. (1), leads to longer ranged site- 
off-diagonal interactions in the array, and a smoother distribution of the externally applied 
voltage drop. Such an extended Hubbard model is no longer integrable via the Bethe ansatz 
technique, but Lanczos direct diagonalization investigations [12| indicate that the physics is 
qualitatively similar to that described here. Disorder, neglected in the present treatment, is 



not found to modify our main conclusions, as confirmed by our Lanczos investigations [12 



The scaling form of the capacitance [Eqs. (8) and (9)] in the Mott insulating phase of the 
system is preserved provided the disorder is not sufficiently strong to lead to a compress- 
ible state, although the correlation length ^ is found to depend on disorder. Similarly, the 
fractional polarization charge plateaus shown in Fig. 3 are robust with respect to disorder, 
though the voltage bias of the steps may be shifted. We also remark that for temperatures 
much larger than the effective charge transfer matrix element, the form of the capaci- 
tance peaks given in Eq. (8) and Fig. 3 will be replaced by a simple derivative of the Fermi 
function, of width fc^T; however, the peak positions can still be used to distinguish between 
metallic and Mott-insulating behavior. 

In conclusion, we have investigated coherent charge transfer in a strongly-correlated 
artificial linear molecule of tunnel-coupled quantum dots. The polarizability in the Mott- 
insulating phase of the system was found to exhibit a universal scaling form analogous to the 
conductivity of the system, while the polarization charge in the dilute metallic phase of the 
system was found to be quantized in units of e/2. Equilibrium charge transfer measurements 
present an intriguing alternative to transport measurements to characterize the electronic 
states of ultrasmall structures. We believe that capacitance measurements carried out in 
carefully fabricated quantum dot chains should be able to observe the novel charge transfer 
resonances and the universal scaling behavior as well as the metallic charge fractionalization 
phenomenon predicted here. 
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FIGURES 



FIG. 1. The quantum corrections to the capacitance plotted in units of /t of a chain of 4 
quantum dots as a function of the bias voltage V . The mimbcr of excess electrons in the chain and 
the interaction strength are as indicated in a legend. Inset: The equivalent circuit of the quantum 
dot array under study. 

FIG. 2. (a) The charge transfer- induced resonant capacitance peaks for Mott-insulating chains 
of 2, 4, and 8 dots (plotted as indicated in a legend), varying by several orders of magnitude in 
height and width, are shown to collapse on the rescaled capacitance peak given by the Eq. (8) 
(solid line), (b) The effective coupling tgjj between the boundary dots of the Mott-insulating 
L-dot chains is plotted for U/t = 6, 7, 8, 9, 10. The solid line with slope minus unity is shown to 
emphasize the scaling of the data given by the Eq. (9). 

FIG. 3. Quantum correction to the polarization charge induced on the external capacitor plates 
versus bias voltage for a Hubbard chain of quantum dots in the metallic phase {N < L). Note 
that Q is quantized in units of e/2 {N <^ L or N ^ L) reflecting a fractional (e/2) charge transfer 
within the chain. 
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TABLES 





I 

2t/e <V <Vi 


II 

Vi<V <V2 


III 

V>V2 


N = L 


kL=iHeV/2t) 
kL-i =Tr-i\n{eV/2t) 


kL=iHeV/2t) 
kh-i = TT + zsin~-^(ism + U /2t) 
Am = sinA:L — iU/4:t 


kL = i\n{eV/2t) 

ki-i = —i shi^^ (i sin k^ + U /2t) 
Xm = sin/cL — iU/At 


N <L 


kN = iln{eV/2t) 


k^r = iln{eV/2t) 
Xm = sin /cjv — iU /At 


kN =iln{eV/2t) 
kN-i = —ism^^{isinkN + U/2t) 
Xm = sin fcjv — iU /At 



TABLE I. Complex roots of Eqs. (4) and (5) (with exponential accuracy as L ^ oo) corre- 
sponding to the bound and antibound states present in the ground state of Eq. (1), as a function of 
the bias V. Here eVi = C//2 + [([7/2)2 +4t2]i/2 and eV2 = U+ [U'^ + 41^]^/"^ . In the insulating phase 
of the model (A^ = L), polarization of the system proceeds via transfer of the antibound state /c^-i 
from one end of the array to the other at the boundary of regions I and II. In the metallic phase, 
polarization of the system proceeds via the successive trapping of electrons on the boundary dot 
with attractive potential at the onset of regions I and III. 
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